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The dispersions, weights, and widths of the peaks of the single-particle spectral function in the 
presence of pair correlations, for a Fermi gas with either attractive or repulsive short-range inter- 
particle interaction, are determined in the normal phase over a wide range of wave vectors, with a 
twofold purpose. The first one is to determine how these dispersions identify both an energy scale 
known as the pseudo-gap near the Fermi wave vector, as well as an additional energy scale related to 
the contact C at large wave vectors. The second one is to differentiate the behaviors of the repulsive 
gas from the attractive one in terms of crossing versus avoided crossing of the dispersions near the 
Fermi wave vector. An analogy will also be drawn between the occurrence of the pseudo-gap physics 
in a Fermi gas subject to pair fluctuations and the persistence of local spin waves in the normal 
phase of magnetic materials. 
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I. INTRODUCTION 

Local order of short-range nature in the normal phase 
of an ultra-cold Fermi gas above the superfluid tempera- 
ture Tc has recently been the subject of intense interest, 
owing to experimental and theoretical advances which 
have hinged on this local order from different perspec- 
tives. 

The experimental interest [H, Q has mostly focused on 
the issue of the pseudo-gap Apg, which is a low-energy 
scale that in these systems evolves in temperature with 
continuity out of the pairing gap present in the broken- 
symmetry (superfluid) phase [3|]. On physical grounds, 
this continuous evolution is due to the persistence of 
"medium-range" pair correlations, which are the remnant 
above Tc of the long-range order below Tc- 

The theoretical interest has been prompted, on the 
other hand, by the introduction of a number of universal 
relations due to Tan 0, [Hj , which are due to the inter- 
particle interaction being of the contact type and affect 
several physical quantities. These universal relations all 
depend on a coupling- and temperature-dependent quan- 
tity named the contact C, which can in turn be con- 
veniently expressed in terms of a high-energy scale Aqo 
0. The fact that C specifies, in particular, the strength 
of "short-range" pair correlations between opposite spins 
implies that, in ultimate analysis, the high-energy scale 
Aoo associated with C and the low-energy scale Apg asso- 
ciated with the pseudo-gap both originate from the same 
kind of pair correlations, which remain active above Tc 
even in the absence of long-range order. 

In this paper, we aim at organizing these two energy 
scales into a single wave- vector dependent function A(fc), 
of which Apg represents the value about the Fermi wave 
vector kp and Aqo its behavior for k much larger than 
kp, corresponding to medium- and short-range pair cor- 
relations, in the order. In practice, from the numerical 
calculations it is meaningful to determine the values of 



A(fc) just in these two intervals, namely, for k ^ kp (ob- 
taining Apg) and k ^ kp (obtaining Aoo)- This wave- 
vector dependence arises even though the inter-particle 
interaction is of the contact type, which at the mean-field 
level below Tc would instead give rise to a wave-vector 
independent gap. 

To this end, we shall analyze in detail the dispersions of 
the peaks of the single-particle spectral function for vari- 
ous couplings across the BCS-BEC crossover and temper- 
atures above Tc, and show how they can rather accurately 
be represented by BCS-like dispersions with a character- 
istic "back-bending" for the occupied states [H, 0] . These 
dispersions will be obtained within the t-matrix approxi- 
mation for an attractive inter-particle interaction follow- 
ing the approach of Ref.Q, which was recently applied 
to account for the experimental data on ultra-cold Fermi 
gases Q. In addition, we will show that the weights of 
the two peaks of the single-particle spectral function can 
also be described by BCS-like expressions. Determining 
these weights will also be useful to obtain the asymptotic 
value of A (A;) for large fc, where tracing the dispersions 
may become ill-defined owing to the strong broadening of 
the large-Zc structure of the single-particle spectral func- 
tion at negative frequencies. 

The importance of determining the weights (besides 
than merely focusing on the existence of the pseudo-gap) 
is in line with the emphasis that was given from the early 
days of the BCS theory of superconductivity to the role 
of the "coherence factors" . Their presence, in fact, made 
the BCS theory soon accepted as the correct one, because 
it was then possible to account for the counter-intuitive 
outcomes of different experiments that could otherwise 
not be understood only on the basis of the occurrence of 
a gap in the single-particle spectrum 0, [l^ . 

The crossed check, between the dispersions and weights 
of the two branches of the single-particle spectral func- 
tion, will therefore represent a fingerprint of the sur- 
vival in the normal phase of typical BCS-like features 
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due to strong pairing fluctuations. Differences, however, 
between the broken-symmetry phase below Tc and the 
pseudo-gap phase above will mostly appear in the 
widths of the peaks of the single-particle spectral func- 
tion, which are much broader above than below as 
expected in the absence of a truly long-range order. 

In this context, we shall also resume an argument that 
was raised in Ref. [lH , according to which the occurrence 
of the above mentioned back-bending for k ^ kp should 
not reflect per se the presence of a pseudo-gap for k about 
kp. This is because the structure at large k, which is 
related to the contact C, can be found even in a normal 
gas with repulsive inter-particle interaction. 

Accordingly, we shall argue that the main differences, 
between the features of the single-particle spectral func- 
tion for a Fermi gas with repulsive or attractive interac- 
tion in the normal phase, appear actually for k about kp. 
Specifically, an avoided crossing results in the dispersions 
of the two peaks of the single-particle spectral function 
in the attractive case, while a crossing occurs in the re- 
pulsive case. In the attractive case, the energy spread of 
the avoided crossing is directly related to the pseudo-gap 
energy scale Apg. On physical grounds, this difference 
between avoided crossing and crossing is a consequence 
of particle-hole mixing, which survives at a local level in 
the attractive case when passing from below to above Tc 
but is absent in the repulsive case. 

The overall purpose here, therefore, is not to establish 
specific criteria for the existence of a pseudo-gap phase in 
a Fermi gas with attractive interaction. Rather, we shall 
be interested in framing the total amount of information, 
which can be extracted from the single-particle spectral 
function of an interacting Fermi gas subject to pairing 
fluctuations above Tc, into a unified picture where analo- 
gies and differences with respect to a simple BCS-likc 
description below Tc can be emphasized. 

It is, nevertheless, relevant to provide at this point an 
(albeit concise) overview of the major relevant work done 
previously by several groups on the issue of the pseudo- 
gap, also to recall how this concept had developed in the 
context of a Fermi gas with attractive interaction. The 
prediction for the existence of a pseudo-gap in the normal 
phase of strongly interacting ultra-cold Fermions was in- 
troduced in Ref. [l2l[ w ithin a two-channel fermion-boson 
model and in Ref.[l3| within a single-channel fermion 
model, before the observation of superfluidity in these 
gases. These works were, in turn, based on earlier studies 
which applied the physics of the BCS-BEC crossover to 
the high-temperature cuprate superconductors. In that 
context, initial work interpreted the normal state of a su- 
perfluid in the crossover regime between BCS and BEC 
as a phase of uncorrelated pairs Q or as a spin-gap 
phase [l3|. Later, it was shown that this phase reflects 
a normal phase pseudo-gap, which displays peculiar fea- 
tures in the fermionic spectral function that reflect the 
presence of a pairing gap in the supcrfluid phase [lit ■ 



Extensive theoretical work on the pseudo-gap issue for a 
Fermi gas with attractive interaction was reported more 
recently in Refs.[l7- 22|. 

Finally, in the present paper a similarity will be high- 
lighted between the pseudo-gap physics resulting from 
pairing fluctuations above Tc and the persistence of spin 
waves over limited spatial regions in the normal phase of 
ferromagnetic (or anti-ferromagnetic) materials. Besides 
being of heuristic value for envisaging the local order as- 
sociated with the pseudo-gap, this analogy evidences how 
the current debate about the occurrence of a pseudo-gap 
in an ultra-cold Fermi gas retraces a similar debate that 
went on for some time about the persistence of spin waves 
in magnetic materials. 

The paper is organized as follows. In Section II, the 
dispersions, weights, and widths of the peaks in the 
single-particle spectral function for the attractive case are 
studied in detail, to determine how the energy scale as- 
sociated with the pseudo-gap about kp evolves for large 
k toward the energy scale associated with the contact 
C. In Section III, the single-particle spectral function 
for the repulsive case is contrasted with that for the 
attractive case, to bring out the issue of the crossing 
vs avoided-crossing of the dispersion relations about kp 
which clearly differentiates between the two cases. In 
Section IV, an analogy is drawn between the pseudo-gap 
physics and the persistence of spin waves in magnetic 
materials, and a suggestion is made for an additional ex- 
perimental evidence for the occurrence of a pseudo-gap. 
The Appendix gives analytic details about the treatment 
of pair fluctuations in the repulsive case to obtain the 
single-particle spectral function over a wide range of k. 



II. THE ATTRACTIVE CASE: PSEUDO-GAP VS 
CONTACT 

In this Section, we consider a homogeneous Fermi gas 
with an attractive interaction VQS{r — r') of short-range 
between opposite spin atoms with equal populations, 
whose strength vq can be eliminated in favor of the scat- 
tering length ap via the relation: 



Airap 



1 



dk m 
(27r)3 k2 



(1) 



Here, m is the particle mass, k a wave vector, and kg a 
wave-vector cutoff which can be let — > oo while vq — > 
in order to keep ap at a desired value (we set h = 1 
throughout). 

Since vq < 0, ap can be positive as well as negative, 
and the dimensionless interaction parameter (kpap)~^ 
ranges from (kpap)~^ ^ —1 in the weak-coupling (BCS) 
regime, to (kpap)^^ > +1 in the strong-coupling (BEC) 
regime, across the unitary limit where \ap\ diverges and 
{kpap)~^ = 0. In practice, the BCS-BEC crossover 
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region of most interest is limited to the interval — 1 ^ 

{kp apy^ < +1. 

In the superfluid phase well below T^, a description of 
the BCS-BEC crossover results already at the mean-field 
level, while in the normal phase above Tc inclusion of 
pairing fluctuations is required to get physically mean- 
ingful results. Pairing fluctuations, in particular, turn 
the characteristic BCS mean-field energy gap below Tc 
into a pseudo-gap above as discussed next. 

A. Mean-field description belov^r Tc 

The simplest description of the BCS-BEC crossover 
results within mean field for temperatures T below 
by supplementing the equation for the BCS gap A 



dk (l-2f{E]^) 



(2^)3 V 2£;k 
with the density equation 

dk 

(2^ 



m 

k2 



All ap 



(2) 



(f-/(i?k))(l-^ 



(3) 



Here, E^, = + with = kV(2m) - ^l and 

f{E) = {e^/C'BT) _,_ jg i-j^g Fermi function {n being 
the fermionic chemical potential and the Boltzmann 
constant). Note that the mean-field gap A does not de- 
pend on fc = |k| owing to the short-range nature of the 
inter-particle interaction. 

When one looks at the structures of the single-particle 
spectral function A(k,uj) within the BCS approximation, 
at a given k two sharp peaks appear centered at the fre- 
quency values 



(4) 



with weights (1 ± ^k/£'k)/2, respectively j23| . 

The dispersion relations ([4]) and the corresponding 
weights are shown in Fig. [1] for three characteristic cou- 
plings across the BCS-BEC crossover at zero temperature 
(dashed-dotted lines). Here, Ep = kp/{2m) is the Fermi 



energy with kp = (Stt n) 
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B. Pairing fluctuations above Tc 

The above picture gets somewhat modified when pair- 
ing fiuctuations beyond mean field arc considered below 
Tc (23 |. It is, however, above Tc that inclusion of pairing 
fluctuations alters mostly the behavior of A{k, ui) from its 
trivial mean-field description with A = 0, whereby only 
a single sharp peak of unit weight survives consistently 
with a Fermi-liquid picture [§] . 

In the context of the BCS-BEC crossover, a non-trivial 
behavior of the spectral function (at and) above Tc results 
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FIG. 1: Dispersion relations (left panels) and corresponding 
weights (right panels) for three different couplings, as ob- 
tained at T = within mean field (dashed-dotted lines) and 
at Tc with the inclusion of pairing fluctuations (diamonds). 
Energies are in units of Ef- 



when including pairing fiuctuations within the t-matrix 
approximation. It is still possible to identify two peaks 
in A{k,Lu) for given k over an extended range of coupling 
and temperature, by locating their positions and deter- 
mining their weights and widths. It is found that the 
positions of these peaks can be rather well represented 
by a BCS-like dispersion of the form (jj]), provided the 
mean-field A is replaced by a "pseudo-gap" value Apg 
that remains finite above Tc- This finding was explicitly 
demonstrated in Ref.0 only for the dispersion relations 
not too far from kp and on the BCS side of the crossover. 

As an example, we report in Fig. [1] the dispersion re- 
lations and weights of the two peaks of A{k, u) at Tc for 
three couplings across the BCS-BEC crossover, obtained 
according to the t-matrix approximation (diamonds) ■ 
In all cases, the similarity with the corresponding values 
obtained for these quantities within mean field at T = 
(dashed-dotted lines) appears striking. 

Marked differences appear instead for the widths of 
the peaks of A{k,LS), when passing from the mean-field 
description below Tc where they are delta-like, to the 
t-matrix description (at and) above Tc where they are 
broad and overlapping. This is shown in Fig. [21 where 
A(k,ui) is plotted versus uj at unitarity for the wave vec- 
tor where the maximum of the lower dispersion relation 
occurs (that is, 0.76fci;' within mean field at T = and 
0.91fci? with the inclusion of pairing fluctuations at Tc). 
This picture also evidences how a real gap at T = trans- 
forms into a pseudo-gap at Tc, through a partial filling 
of the spectral function in the region between the two 
peaks. In the present paper we shall dwell extensively on 
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FIG. 2: Single-particle spectral function A{k,uj) versus lj at 
unitarity, within mean field at T = (dashed-dotted lines) 
and with the inclusion of pairing fluctuations at Tc (full line), 
for the wave vector where the maximum of the lower disper- 
sion relation occurs in both cases (see text). 



this and related ideas. 

The t-matrix approximation that we adopt in this pa- 
per to obtain A{k, uj) above Tc corresponds to the follow- 
ing choice of the fermionic self-energy Q: 

S(k,w„)= - y^fcfiT^ ro(q,l}.) 

X Go(q- k,0,^ - a;„) (5) 

where u;„ = {2n + VjuksT (n integer) and fli, = 2i>T:kBT 
[v integer) are fermionic and bosonic Matsubara frequen- 
cies, in the order, Go(k, a;„) = (ia;„ — ^k)~^ is the bare 
fermionic single-particle Green's function, and Tq{(\^VL^) 
is the particle-particle ladder given by 



ro(q, a. 



(2^ 



ksT'^ Go(k, iu„) 



(6) 



X Go(q-k, ^ k2 



The single-particle spectral function is then obtained 
through analytic continuation iun lu + irj to the real 
frequency axis (77 = 0+): 



1 



ImE(k,a;) 



- - RcS(k, uj)Y + [Iml](k, uj) 



(7) 

The shape of A{k, uj) versus uj thus depends crucially on 
the interplay between ReS(k, cj) and Iml](k, w) for the 
chosen value of k = |k|. 

A derived quantity of interest is the single-particle den- 
sity of states: 



N{uj) 



Aik,uj) 



(8) 



The averaging that this definition introduces on A(k, uj) 
over an extended range of k can be of support to the 



presence of a pseudo-gap, in cases when the two peaks of 
A(k, uj) strongly overlap just in the range of k where the 
two branches of the dispersion come close to each other 
(cf. Fig. [T|). In these cases, in fact, a strict definition 
of the pseudo-gap as a depression of the spectral weight 
just in this range of k would lead one to conclude that 
pseudo-gap phenomena were absent in the single-particle 
excitations, while they still appear clearly over a more 
extended range of k. 

C. Inputs from experiments on ultra-cold Fermi 
atoms 

The original motivation for looking at A(k, uj) has 
been the issue of "preformed pairs" in high-temperature 
(cuprate) superconductors, before the occurrence of the 
BCS-BEC crossover was explicitly demonstrated with 
ultra-cold Fermi atoms (cf., e.g., Ref.^25i]). 

In this context, the interest in the detailed shape of 
A{k, ui) above Tc across the BCS-BEC crossover has con- 
siderably raised lately, after a new measurement tech- 
nique was introduced to probe directly the single-particle 
excitations of a Fermi gas Intensity maps were thus 
obtained for the single-particle excitation spectra, relat- 
ing the single-particle energy to the wave vector. More 
recently, new measurements performed over an extended 
temperature range above Tc jj] have revealed a BCS-like 
dispersion with a characteristic "back-bending" close to 
kp which identifies a pseudo-gap energy scale and per- 
sists well above Tc. 

This finding gives us motivations for extending the the- 
oretical analysis of Ref . 0] for the dispersions of the peaks 
of A{k,uj) across the unitary region, although the widths 
of the peaks can increase considerably with respect to 
the BCS side, reflecting the fact that quasi-particle exci- 
tations may be poorly defined. Specifically, the combined 
experimental and theoretical analysis of Ref. Q suggests 
us to concentrate our efforts in the coupling range ap- 
proximately between [kpap)"^ = and {kpap)~^ = 0.4. 

D. Emergence of the contact in A{k,uj) 

Yet, it was pointed out [llj that the persistence of 
the back-bending for large k (^ kp) is dominated by 
interaction effects that do not reflect the pseudo-gap close 
to kp. Rather, it is connected with the universal tail 
of the wave- vector distribution n(k) of a dilute Fermi gas, 
whose coefficient is given by the Tan's contact G 0, 

This property can be readily verified within the t- 
matrix approximation that we use to obtain A{k,uj). 
When k'^/{2m) or |cl>„| are much larger than the energy 
scales ksT and |ju|, in fact, the self-energy (O can be 
approximated by: 



S(k,a;„) 



;n/ro(k,w„) - A^Go(k,-w„) . (9) 
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Here, 



with 



(2^ 



(10) 



nf{k)^kBTY,e"^"''Go{Ki^n) (11) 

n 

is the free density associated with Co for given ^, and 

is the square of the high-energy scale introduced in Ref. Q 
that was mentioned in the Introduction. The two terms 
on the right-hand side of the approximate expression ([9]) 
originate from the singularities in the complex frequency 
plane of the single-particle Green's function Go and of the 
particle-particle ladder Tq, in the order, once the sum 
over the Matsubara frequency in the expression ([5]) of 
the fermionic self-energy is transformed into a contour 
integral. 

Analytic continuation za;„ — )■ to + ir] to the real fre- 
quency axis then results into the following approximate 
expression for large k: 



A? 



A? 



Aik, 0.) ^ (^1 - ^ j 5(0. - Ck) + + (13) 

which presents indeed a well-defined structure at the neg- 
ative frequency u! = — ^k- One obtains correspondingly: 



n{k) 



-t-CxD 



did f{Lo) A{k,uj) 



(m Ac 



Hi 



fc4 



(14) 



yielding the relation C = (mAoo)^ between the contact 
C and Aoo, which will be extensively used below. 

In practice, the structure of A{k,uj) for negative real lj 
at large k is not delta-like but spreads over a sizable fre- 
quency range. This difference from the approximate re- 
sult (jl3p stems from the non-commutativity between tak- 
ing the analytic continuation and performing the large-A; 
expansion of the self-energy, as recalled in the Appendix. 
Nevertheless, the actual structure oi A{k,u!) for negative 
u! preserves the same total area A^/(4^^) found above 
in the expression ^ 



E. Connecting the two energies Apg and Aoo 

With these premises, it seems natural to frame the 
low-energy scale Apg and the high-energy scale Aoo into 
a unified physical picture, in which they emerge from 
A(k, uj) in the two distinct ranges of wave vectors k kp 
and k ^ kp, respectively. To this end, we shall extend 
to the unitary limit and beyond the analysis that was 
limited in Ref.Q to the BCS side of the unitary region. 



by following the dispersions, weights, and widths of the 
peaks of A{k^uj) from = 0, through k « kp, and up 
to A: 3> kp, even in cases when these peaks appear quite 
broad and overlapping. 

On physical grounds, the evolution, from Apg when 
k ~ kp to Aoo when k ^ kp, is expected on the basis 
of a (local in space and transient in time) order which is 
established by pair fluctuations above Tc, in the absence 
of long-range order (as described by mean field below 
Tc). This is in line with the definition of the contact 
C (and thus of Aoo) through the short-range behavior 
of the pair-correlation function between opposite spins 
[1], which corresponds to k ^ kp; while the pseudo-gap 
Apg is expected to depend on pair-correlations that are 
established more extensively over medium range, which 
corresponds to fc « fc^- 

As a consequence, we expect the "pseudo-gap physics" 
to be associated with pair correlations which are built 
at intermediate distances of the order of kp^, while the 
"contact physics" with pair correlations that survive even 
at smaller distances (<C kp^). Both quantities Apg and 
Aoo are thus affected by the same sort of pair corre- 
lations in the particle-particle channel, to which the t- 
matrix approximation that we adopt in this paper pro- 
vides an important contribution. This is because the 
particle-particle ladder propagator [given by Eq.® in 
the attractive case and by Eq. ip^ in the repulsive case] 
generalizes to a many-body environment the two-body t- 
matrix, which describes two-body binding and scattering 
of unbound particles at the same time. 

Within this approximation, (the square of) Aoo is de- 
fined from Eq. (|12p as a wave- vector and frequency av- 
eraging of the particle-particle ladder propagator. The 
wave-vector and frequency structures of the same prop- 
agator give also rise to the pseudo-gap Apg, which then 
emerges as a characteristic low-energy scale in the single- 
particle excitations. 

Consistently with this physical picture of locally estab- 
lished pair correlations, we expect Apg to survive above 
Tc over a more limited temperature range than Aoo, since 
thermal fiuctuations act first to destroy the order estab- 
lished over intermediate distances. Our analysis about 
the characteristic features of the spectral functions will 
accordingly be extended over a meaningful temperature 
interval above Tc- 

F. Working procedures 

Working experience on the single-particle spectral 
function suggests us to identify the low-energy scale Apg 
in the range ^ fc ^.kp, while the high-energy scale Aoo 
can be extracted with sufficient accuracy already from 
the not too extreme range 2kp ^ fc ^ 4kp. In a more 
extreme range of k (> 'ikp), in fact, it would become 
quite difficult to determine A{k,uj) for large negative uj. 

Owing to the shape of the w-structures of A{k,uj) for 
given k, different strategies need to be adopted in the 



6 



above two ranges of wave vectors. Namely, in tlie range 
^ fc < 2fc_F the sfiape of the two peaks of A{k,u}) per- 
mits, in practice, to both follow their dispersions for vary- 
ing k and identify their weights in terms of the total area 
they comprise. In the range 2kp ^ k 4^^, on the 
other hand, the structure of A(k, lo) at negative frequen- 
cies is so broad that only its total area can be reasonably 
identified. 

Range < k < 2kF: 

Let's first consider the range < fc < 2fci?. Here, the 
dispersions that wc are able to determine independently 
for the two peaks at positive and negative frequencies are 
fitted, respectively, by the BCS-like expressions: 



W(±)(fc) = ± 




pg(±) 



(15) 



where a different pseudo-gap energy Apg is introduced for 
the upper (+) and lower (— ) branches. The fitting identi- 
fies, in addition, the locations k^i^^) for the "up-bending" 
of the upper branch and fcL(-) for the "down-bending" 
of the lower branch [2^. The fittings are carried out by 
a x^-analysis. We have found that fcL(+) is consistently 
smaller than fcL(-) in all situations we have examined. 

Because the peaks at negative and positive frequen- 
cies are in general broad and partially overlapping with 
each other, we have chosen to determine operatively their 
weights at a given k by integrating A{k^Lo) from — oo up 
to a; = and from w = up to -l-oo, in the order. We 
then fit these values obtained for several k by the BCS- 
like expressions: 



v{kf = \\ l- 



2^ ~ A^cff 



^2 

Pg(-)' 



for the lower branch at negative frequencies, and 



(16) 



u{kf = i I 1 + 



, 2m 



Moff ) 



A2 

pg(+) ' 



(17) 



for the upper branch at positive frequencies. Here, 
^pg(±) are the pseudo-gap energies already determined 
from the fitting (jlSp to the dispersions, and /^cff is a new 
parameter common to the two branches which is deter- 
mined from the position where the weights cross. We 
shall find that the value of y/2infics is intermediate be- 
tween and fcL(-) obtained from the dispersions ([15]) 
of the two branches, consistently with enforcing a com- 
mon value fj,es in Eqs. ([T6| and ([TTl). 

Note that the expressions (fTH]) and (fTT]) do not require 
a priori the sum u{k)'^ + v{k)'^ to be unity, as it would be 
the case for the coherence factors in a strict BCS descrip- 
tion. Deviations of the sum of the expressions (jl6p and 



pTl) from unity can thus be taken as a test for the validity 
of the effective BCS description that we are attempting 
to establish in the normal phase close to T^. 

Concerning, finally, the widths of the peaks of ^(fc,a;), 
they will be conventionally determined as the full widths 
at half maximum. Whenever necessary, however, these 
values will be compared with the numerical values ob- 
tained alternatively by a two-Lorentzian fit to the peaks 
of A{k,ijj), a procedure which is of course more reliable 
the more these peaks are separated in frequency. 

The BCS-like expressions (US]), dH]), and ^ are 
meant to test the persistence of a BCS-like description 
in the normal phase of a Fermi gas with an attractive 
pairing interaction. In addition, the presence in Eq.()15|) 
of the "Luttinger wave vectors" fci(±) highlights the per- 
sistence of an underlying Fermi surface for the single- 
particle excitations, which represents the last remnant of 
what would be a Fermi-liquid description of the Fermi gas 
if attractivepairing interactions above Tc would not be 
considered [8[. Nevertheless, the large widths associated 
with the peaks of A{k,u;) represent per se an evidence 
that a Fermi-liquid description above Tc does not apply 
in the presence of pairing fluctuations. 

Range 2kF < k < 4kF: 

In this range, the structure at negative frequencies in 
A{k, uj) becomes so spread and broad that it is meaning- 
less to determine its dispersion numerically and then try 
to fit it by an expression similar to ()15|) . In this case, 
however, it remains meaningful to determine the total 
area of the broad structure at negative cj over a chosen 
mesh of k and then make a x^"fit to these values through 
the following expression which is inspired by Eq. (|13p : 
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(18) 



where Aiargc and /^largc are fitting parameters to be de- 
termined in this range of "large" k. 

In practice, it is convenient to set /xiarge at the corre- 
sponding value of the thermodynamic chemical potential 
fi from the outset (thus leaving Aiarge as the only fit- 
ting parameter). This is because in the coupling range of 
interest ^large is small enough that it becomes meaning- 
less to extract it from the denominator of Eq. (|T8)) where 
k'^ /(2m) dominates in this range of k. 



We will check whether the value of Ai 



determined 



in this way coincides with the value of Aoo obtained in- 
dependently by the expression (fT^ . and how it differs 
from the value Apg(_) of the pseudo-gap obtained above 
near the region of the back-bending of the lower branch. 

G. Results for dispersions, weights, and widths 

We pass to determine the quantities of interest accord- 
ing to the procedures outlined above. We shall specif- 
ically consider the two coupling values (kFaF)~^ = 
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FIG. 3: Dispersions (upper panel) and corresponding weights 
(lower panel) at unitarity and T — T^. Circles (squares) and 
full (dashed) lines represent the results of the numerical cal- 
culation and of the BCS-like fits for the lower (upper) branch 
at negative and positive frequencies, respectively. Energies 
are in units oi Ep- 



and {kFa,F)~^ = 0.25 as representatives of the coupling 
range where pseudo-gap phenomena are expected to be 
maximal. Two representative temperatures will also be 
considered for each coupling. 

Figure [3] shows the dispersion relations and weights of 
the two peaks of A{k,u)) at Tc when (kpap)^^ = in 
the range Q < k < 2k p, as obtained from the numerical 
calculation based on Eq. ([7]) and from the fits obtained ac- 
cording to Eqs. (fT5l) - ()17p . In this case the fitting param- 
eters are fci(_) ~ 0.78fcF and Apg(_) = 0.83Ef for the 
lower branch, and fcL(-(-) = 0.62fci? and Apg(-^_) = O.TAEp 
for the upper branch, while fXoS = OAlEp is quite close 
to the corresponding value of the thermodynamic poten- 
tial ^ = 0.365Ep obtained within the t-matrix approxi- 
mation. [The results of the numerical calculations have 
already been reported in the central panels of Fig. [1] 
although with the different purpose of comparing them 
with the mean-field description at T = 0.] 

In this case the BCS-like fits are excellent for the dis- 
persions and quite good for the weights, especially near 
the value k = O.GAkp where the weights exchange with 
one another. Note also that, while the numerical values 
of the weights for each value of k are specular to each 
other about one half, the fitted values are not always so 
indicating deviations from their sum being unity. 

This success of a BCS-like interpretation for the disper- 
sions and weights of the peaks should be complemented 
by the further information about their widths. This is 
done in Fig. 21 where in the upper panel the shape of 
A{k,uj) at Tc and (kpap)^^ = is shown explicitly for 



FIG. 4: Unitarity limit and T = Tc. Upper panel: A{k,w) 
vs uj for the wave vectors k = (0.6, 0.7, 0.8, 0.9, 1.0, 1.1, 1.2)kF 
corresponding to the peaks at negative frequency from top 
to bottom (here full and dashed lines alternate to help the 
analysis of the figure). Lower panel: Widths (in units of 
Ef) of the peaks at negative (circles) and positive (squares) 
frequencies. 



several wave vectors, while in the lower panel the cor- 
responding widths of the peaks at negative and positive 
frequencies are reported over a wider set of wave vectors. 

In all cases, the widths are rather large (being com- 
parable to Ep) and show strong deviations from what 
would be expected for a Fermi liquid picture, according 
to which they should acquire a minimum value at about 
kp. These deviations from a Fermi liquid picture are 
of course expected for a Fermi gas with attractive inter- 
action, taking further into account that at unitarity the 
value of Tc is a considerable fraction of the Fermi temper- 
ature Tp whereas a Fermi-liquid description holds only 
for T <^Tp^. 

The above analysis of pseudo-gap phenomena around 
kp is expected to remain meaningful for temperatures 
larger than Tc, but not exceeding the pair-breaking tem- 
perature scale T* where a "preformed-pair scenario" is 
bound to fade away. In particular, at unitarity the value 
of T* (as estimated by the mean-field critical temper- 
ature) is about twice the value of Tc given by the t- 
matrix approximation we are considering [isl ] [l3| . As 
a representative case of a temperature above Tc, Fig. [5] 
shows the dispersions and weights obtained from A{k, oj) 
at unitarity and T = 1.2 Tc- The fitting parameters are 
now = 0.99fcF and Apg(_) = QASEp for the lower 



branch, and k 



Li+) 



Omkp and A 



pg(+) 



= 0.62Ep for 



the upper branch, while /icff ~ OASEp (to be compared 
with the thermodynamic value fi — 0.39Ep). 



FIG. 5: Dispersions (upper panel) and weights (lower panel) 
at unitarity and T = 1.2 Tc. Conventions are as in Fig. [31 



Note that in this case the analysis of the dispersion 
of the lower braneh had to be interrupted over a non- 
negligible interval of k about kp, because in this interval 
the structure of A{k, uj) at negative frequencies is almost 
completely masked by the stronger structure at positive 
frequencies. This represents a signal that pseudo-gap 
phenomena are beginning to fade away at this temper- 
ature. We have nevertheless performed a BCS-like fit 
to the part of the dispersion that can still be clearly 
identified, as shown by the full curve in the upper panel 
of Fig. [5l By our procedure no problem instead arises 
in identifying the corresponding weights reported in the 
lower panel of Fig. [5l which follow again a BCS-like dis- 
persion although with less accuracy than those shown at 

in Fig. El 

The corresponding shapes of A{k^uj) vs ui for a cho- 
sen set of k across kp are shown explicitly in the upper 
panel of Fig. [6l from which one can appreciate the phe- 
nomenon mentioned above, when the structure oi A{k, ui) 
for the lower branch becomes a shoulder attached to the 
structure of the upper branch. The corresponding broad- 
enings of these two structures are reported in the lower 
panel of Fig. [6l which reinforces our conclusion about the 
non-Fermi-liquid nature of the system. 

A question naturally arises, about whether or not these 
profiles of A{k,ijj) still allow one to identify the presence 
of a pseudo-gap in the crucial range of wave vectors about 
kp. As the upper panel of Fig. [5] shows, in fact, an overall 
BCS-like fit to the lower branch can be attempted even 
in this case, because the two structures of A{k, uj) remain 
distinct from each other away from kp. 

The relevance of this restricted interval about kp can 
be strongly reduced by averaging the profiles of A{k, to) 
over all wave vectors, in the way it is done in the def- 
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FIG. 6: Unitarity limit and T = 1.2 T^. Upper panel: A{k,uj) 
vs UJ for the same wave vectors as in the upper panel of Fig. 2] 
Lower panel: Widths (in units of Ef) of the peaks at negative 
(circles) and positive (squares) frequencies. 
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FIG. 7: Density of states per spin component vs u calculated 
at unitarity within the t-matrix approximation for the tem- 
peratures: T ^ Tc (fuU line), T = 1.2 Tc (long-dashed line), 
T = 1.4Tc (short-dashed line), and T = 1.65 Tc (dot-dashed 
line). In the present case, T* 2Tc. The dotted line shows 
the corresponding mean-field result when T — 0. The non- 
interacting value mkp /{2-k'^) of N{uj = 0) at T = is used 
for normalization. 



inition ([5]) of the single-particle density of states N{uj). 
Figure [7] shows a plot of A''(a;) vs w at unitarity for sev- 
eral temperatures at and above Tc [1^. For increasing 
T, the depression of N(ui) near a; = well survives for 
T ~ 1.2 Tc at which the dispersion of the lower branch 
near kp in Fig.[5]had to be interrupted, and progressively 
disappears for temperatures somewhat below the pair- 
breaking temperature scale T*. That the depression of 
density of states survives at temperatures higher than the 
crossover temperature where the pseudo-gap features dis- 
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FIG. 8: Dispersions (upper panel) and weights (lower panel) 
at T = Tc for the coupling (/cfAf)^^ = 0.25. Conventions are 
as in Fig. [3] 



appear in the spectral function was previously discussed 
in Refs.[ll[lil. At about T*, N{lu) for w = coincides 
(within a few percent) with its non-interacting value eval- 
uated at the same temperature and chemical potential, 
indicating that all effects of pairing have faded away at 
w = (although they will persist at higher temperatures 
for uj <C —Ep, indicating the survival of the "contact" 
even at quite high temperatures jsO])- The density of 
states obtained within mean field at zero temperature is 
also reported for comparison in Fig. [71 and shows two 
sharp peaks located at ±A with A = 0.69 £'f- 

It is important to extend the above analysis past the 
unitarity limit to the BEC side of the crossover (but still 
before the pseudo-gap turns into a real gap associated 
with the binding energy of the composite bosons which 
form in the BEC limit). To this end, Fig.|S]shows the dis- 
persions and weights at Tc for the coupling {kpap)~^ = 
0.25, together with the corresponding BCS-likc fits where 
now fcL(-) = 0.77kp and Apgj^) = 1.09Ep for the lower 
branch, and fcL(+) = 0-28kp and Apg(^-) ~ 0.91Ep for 
the upper branch, while = 0.09Ep (which in this 
case almost coincides with the thermodynamic value). 
Compared with Fig. |3l the dispersions have now become 
quite flat in the range < k kp, while the two weights 
cross at smaller value of fc. At even stronger couplings, 
the lower dispersion down-bends and the upper disper- 
sion up-bends already at fc = while the weights always 
remain well separated from each other for all k (as it 
is shown in the lower panels of Fig. [T] for the coupling 
(kpap)-'^ = 1.0). 

For completeness. Fig. |9] shows the corresponding 
shapes of A{k^uj) across kp (upper panel) as well as the 
broadcnings of two structures of A{k,uj) (lower panel). 



FIG. 9: Coupling (kpaF)''^ = 0.25 and T = T^. Upper panel; 
A{k,LL!) vs oj for the same wave vectors as in the upper panel 
of Fig. [4] Lower panel: Widths (in units of Ef) of the peaks 
at negative (circles) and positive (squares) frequencies. 




FIG. 10: Dispersions (upper panel) and weights (lower panel) 
at T = lATc for the coupling {kpap)^^ — 0.25. Conventions 
are as in Fig. [S] 



When the coupling increases toward the BEC 
regime, the pair-brea king temperature T* increases more 
markedly than Tc [28 1 13 1 and pairing fluctuations are ac- 
cordingly expected to affect A{k, uj) over a progressively 
wider temperature range above Tc- We then report in 
Fig. [TU] the dispersions and widths of the two branches 
of A{k,LL!) for the coupling {kpap)~^ = 0.25 and the 
higher temperature T = 1.4 Tc. Again, a signal that the 
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FIG. 11: Coupling {kpaF)'^ = 0.25 and T = 1.4 Upper 
panel: A{k, lj) vs lj for the same wave vectors as in the upper 
panel of Fig. [l] Lower panel: Widths (in units of Ef) of the 
peaks at negative (circles) and positive (squares) frequencies. 



pseudo-gap is beginning to fade away emerges from the 
analysis of the dispersion for the lower branch, which has 
to be interrupted about kp. The fitting parameters are 
now ki^(^^) = LOSfci? and Apg(_) = Q.QOEp for the lower 
branch, and fcL(+) = Q.2bkp and Apg(+) = Q.^QEp for the 
upper branch, while /Zeff = Q.QQEp (to be compared with 
the thermodynamic value ^ = Q.lAEp). The correspond- 
ing shapes of A(fc,a;) across kp and the broadenings of 
two structures of A{k, ui) are shown, respectively, in the 
upper and lower panels of Fig. [11] 

Beginning with Fig. [21 we have often emphasized that 
one of the major characteristics of the two structures of 
A(fc, uj) (at and) above Tc is their substantial broadening, 
which may hinder in practice a straightforward identifi- 
cation of the pseudo-gap about kp in cases when these 
structures strongly overlap with each other. In these 
cases, however, one may resort to a two-Lorentzian fit 
of the two structures of A{k, uj) which helps separating 
them. This is shown Fig. [T^j where the dashed lines rep- 
resent the two Lorentzians. For instance, by this type 
of fit our previous estimates for the weight (0.29) and 
width {OMEp) of the structure of A{k,uj) at fc = O.Ofcf 
corresponding to the lower branch are replaced by 0.28 
and O.llEp, in the order. Comparable deviations are ob- 
tained in the other cases. These results thus confirm the 
validity of our previous analysis where the weights and 
widths were extracted from A{k, uj) in a simpler fashion. 

The numerical values of the dispersions, weights, and 
widths that were reported in the previous figures were all 
obtained from the detailed profiles of the single-particle 
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FIG. 12: Two-Lorentzian fits (dashed lines) of A{k,uS) vs 
UJ (full lines) at unitarity and T = Tc, when k/kp = 
(0.6, 0.9, 1.2) from top to bottom. 



spectral function A{k,uj), which were also shown in the 
same figures. It may also be of use, however, to organize 
the spectra of A{k,uj) for a range of k and uj into a single 
intensity plot. This is done in Fig.fTHla) for the same set 
of temperatures and couplings considered in the previous 
figures. Similar intensity plots were presented in Refs.|18| 
and . Note that the log scale, used here like in the 
experimental works [H, Q , makes the back-bending more 
evident when compared with the intensity plots presented 
in Refs. 18| and [19[. For the sake of comparison with 
those references, we also report in Fig. [T3lb) the same 
intensity plots in a linear scale. 

Thus far we have concentrated our attention to the 
range ^ k ^ 2kp where the pseudo-gap physics man- 
ifests itself. We pass now to discuss the more asymp- 
totic range 2kp < k < Akp where the contact physics 
emerges. To this end, we adopt the procedure outlined 
in sub-section II-F and determine the parameter Aiargo 



from the expression (|18|) with fiurgc fixed at the corre- 
sponding value of the chemical potential. 

Figure [HI shows the weights of the structure of A{k, uj) 
at negative uj for the values of coupling and tempera- 
ture considered so far, as determined numerically (circles) 
over a mesh of values of k in the range 2kp < k < 4kp 
and then fitted (full lines) in terms of the expression (|18p . 

These fits are also compared with an expression of the 
form (jl8p . where now the low-energy scale Apg(_) that 
was previously determined in the range ^ fc ^ 2fcj? re- 
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FIG. 13: Intensity plots for the single-particle spectral func- 
tion at given temperature and coupling, in (a) logarithmic and 
(b) linear scale. The thin white line identifies the dispersion 
of the lower branch as determined in the previous figures. 



places Aiargc (dashed lines). The appreciable deviations 
from the numerical values of the weights that result show 
that the high-energy scale Aqo can be distinguished from 
the low-energy scale Apg(_) by inspecting the shape of 
A{k,uj) in different ranges of k. Note in particular that, 
as soon the temperature is increased above Tc, Apg(_) 
becomes rapidly smaller than Aoo. This is consistent 
with our expectation that Aqo , being associated with lo- 
cal pair correlations of shorter range with respect to Apg, 
survives at higher temperatures. 

A direct comparison of the temperature dependence of 



Aiargc and Apgj-.-) is shown in Fig. [15] for the two cou- 
plings previously considered. Here, squares represent the 
values of Aiarge obtained from the fittings reported in 
Fig. [m circles are the values of Aoo obtained indepen- 
dently from the expression (jl2[) . and triangles are the 
values of Apg(_) determined from the fittings to the 
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FIG. 14: Weights of A{k, lu) at negative uj (circles) are fitted 
according to the expression discussed in sub-section II-F (full 
lines) for different couplings and temperatures. Dashed lines 
correspond to what would be obtained by using in that ex- 
pression the numerical values of the low-energy scale Apg(_). 
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FIG. 15: Temperature dependence of Aiargc (squares), Aoo 
(circles), and Apg(_) (triangles) in units of Ep for the cou- 
plings (fcFfflF)^^ = (upper panel) and (/cfUf)"^ = 0.25 
(lower panel). 



dispersions. 

It is evident from this figure that Apg(_ 



is a much 



faster decreasing function of temperature than Aqq, 
which reflects the slow decay of the contact C at high 
temperature [s^ . Note also that some discrepancy arises 
between the values of Aiarge and Aoo at increasing tem- 
perature. This is due to the fact that for increasing 
temperature the interval of k from which Aoo can confi- 
dently be extracted should be centered progressively at 



12 



a larger value of k, while in Fig. [15] we have kept it at 
2kp < k < ikp for all temperatures. In any case, the 
difference between Aiargo and Aqo is significantly smaller 
than that between and Apg(_). Note finally that at 
unitarity Apg(_) and Aqo almost coincide with each other 
close to Tc- In this case, a single value can be effectively 
associated with the two energy scales. 

However, the two energy scales soon deviate from each 
other not only for increasing temperature above T^., but 
also away from unitarity on the two sides of the crossover. 
For instance, in the BCS regime A^o = 2Tr\aF\n/'m for 
T < (ma^p)^^ while the pseudo-gap would be exponen- 
tially small in the coupling parameter (/ci^Oi?)"^- In the 
BEC regime, on the other hand, Aqo = \/ ^im/ {rn^aF) 
for T < [ma^p)^^ while in this case the "real" gap in the 
single-particle excitations would equal half the value of 
the binding energy (ma^p)~^ of a composite boson Q. 

It should be mentioned in this context that, by the 
alternative t-matrix approach of Ref. [l^ , a trace of the 
pair-fiuctuation propagator (quite similar to Ea.p2|) for 
A^) was interpreted as representing (the square of) a 
pseudo-gap energy for all couplings and temperatures 
above Tc, thus making in practice the high-energy scale 
Aoo and the low-energy pseudo-gap Apg to coincide with 
each other. This marks a difference between the approach 
of Ref. [l^ and the present one, which keeps instead the 
two energy scales Aoo and Apg distinct from each other. 

III. THE REPULSIVE CASE: CROSSING VS 
AVOIDED CROSSING 

We pass now to consider the occurrence of the two 
energy scales Apg and Aoo from a different perspective, 
which emphasizes the differences one finds near kp for the 
two branches W(-|-)(fc) when considering a Fermi gas with 
attractive or repulsive inter-particle interaction. These 
differences are related to the presence near kp of a finite 
or vanishing value of Apg, while the behavior of uj(^±^{k) 
for k ^ kp which is related to Aoo remains essentially 
the same in the two cases. 

Accordingly, we shall contrast the behavior near kp 
that will result from the single-particle spectral function 
A(fc, uj) for the repulsive case, with that identified already 
in the previous Section for the attractive case. In this 
way, we shall significantly extend the discussion given in 
Ref.fm on A{k,uj) for the repulsive case, that was there 
considered only for large values k ^ kp. 

A. Specular comparisons 

To make a meaningful comparison, we consider the at- 
tractive case with ap < and the repulsive case with 
ap > for the same value {kp\ap\)^^ of the dimension- 
less coupling. In the Appendix we provide the necessary 
analytic details for the less familiar repulsive case, at the 



level of the t-matrix approximation that we use in the 
numerical calculations. 

As we have already discussed, one can not only de- 
termine the dispersions of the two peaks of A(k,uj) but 
also keep track of their weights. In the attractive case 
considered in Section II, this combined information has 
resulted in the phenomenon of avoided crossing which 
typically occurs when two quantum levels with the same 
symmetry evolve as a function of a parameter (sij. In 
this case, pairing fiuctuations induce above Tc the same 
kind of particle-hole mixing which is characteristic of the 
BCS theory below Tc [§]. This mixing, in turn, makes the 
two branches of A{k, lu) to share the same symmetry (be- 
ing partially particle- like and partially hole-like) , in such 
a way that no crossing of the two dispersions occurs (the 
region of their minimum approach being, by definition, 
associated with the pseudo-gap). 

The phenomenon of avoided crossing has to be con- 
trasted with what happens instead in a Fermi liquid, a 
system where particle and hole excitations do not mix. 
In this case, the two branches of A(k,io) are expected 
to cross each other at kp , where they abruptly exchange 
their weights in a similar fashion to what occurs typically 
for a crossing jsij . 

As emphasized in Ref.[ll|, when k ^ kp the occur- 
rence of two branches in A{k,ui) in the place of a single 
one even for a Fermi liquid stems from the requirement 
that for k ^ kp the wave- vector distribution n{k) at zero 
temperature has a tail oc C/k'^, in accordance with Tan's 
argument 0, 0- However, the branch at negative uj has 
an extremely small (albeit non-vanishing) weight, as we 
shall explicitly verify for a Fermi gas with a short-range 
repulsion. 

B. Working procedures 

It is discussed in the Appendix that an appropriate 
choice of the parameters k^/kp, {kpap)~^, and mkpVQ 
entering Eq.([T]) has to be made for the repulsive case. 

In particular, when exploring the region k ^ kp in 
order to extract the quantity Aoo , the values of k should 
not be smaller than, say, Akp (a value consistent with 
the plots reported in Fig. 3 of the second of Refs.[30{. 
where the contact C was calculated numerically within 
the t-matrix approximation). If we choose mkpvo not 
larger than 10 for speeding up the summations over the 
Matsubara frequencies, from Eq. p^ of the Appendix we 
obtain (kpap)^^ ~ 4.4 to be an optimal value of the 
coupling to the purpose. These values for the param- 
eters have also been used in Fig. [^T] of the Appendix, 
where the area of the peak of A{k,u!) at negative en- 
ergies yields the value Aoc/ Ep = 0.1102. As expected 
in this rather extreme weak-coupling regime, a corre- 
sponding calculation made for the attractive case with 
coupling (fci?ai?)~^ = —4.4 yields the comparable value 
Aoo/£^F = 0.0860 [cf. Eq.©]. 
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On the other hand, when exploring the region k k, kp 
to focus on the issue raised above about "crossing vs 
avoided crossing", in the repulsive case the coupling 
{kpap)~^ should not exceed, say, the value 1.5. Other- 
wise, in the corresponding calculation for the attractive 
case it would be extremely difficult to detect numerically 
the occurrence of an avoided crossing when {kpap)^^ 
becomes smaller than —1.5. When (kpap)^^ = 1-5, 
the value of mkpVQ has to be pushed up to 20 to get 
kvi/kp{= 1.37) larger than kp as required. 

C. Results for dispersions, weights, and widths 

Figure [TBI compares dispersions, weights, and widths of 
the two structures of A(fc,w) close to kp, as obtained in 
the repulsive case with {kpap)~^ = 1.5 and T = (left 
panels) and in the attractive case with (kpap)^^ ~ —1.5 
and T = Tc = Q.OblQTp (right panels). This compari- 
son highlights and contrasts the essential characteristics 
found in A{k, to) for a Fermi liquid with no pseudo-gap 
(left panels) and for a non- Fermi liquid with a pseudo-gap 
(right panels), at the corresponding values of {kp\ap\)'^ . 
Note in particular that: 
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FIG. 16: Dispersions, weights, and widths extracted from 
the two features of A(k,u}) near kp for a Fermi system with: 
{kpaF)"^ = 1.5, mkFVo = 20, fco = 1.37fcF, and T = (left 
panels); {kpapy^ = -1.5 and T ^ Tc = 0.0576Tf (right 
panels). Energies are in units of Ep. Full (dashed) lines 
correspond to the structure of A{k,ijj) aX uj < (lj > 0) . The 
meaning of the double arrow is explained in the text. 

(i) The crossing at fc^? of the dispersions of the two struc- 
tures of A{k, uj) in the repulsive case (left-upper panel) 
contrasts with the avoided crossing in the attractive case 



(right- upper panel), whereby the two branches exchange 
their role and remain separated by the amount 2Apg (in 
the present case, the wave vector fci at which the avoided 
crossing occurs between the two branches in the attrac- 
tive case is quite close to kp); 

(ii) The behavior of the spectral weights associated with 
the two structures of A{k,uj) (for ui < and uj > 0, 
respectively) shows an abrupt exchange at kp in the 
repulsive case (left-middle panel) which is typical of a 
level crossing, while a smooth evolution over a spread 
Sk such that (5fc^/(2m) « Apg results in the attrac- 
tive case (the size of Sk is represented by a double ar- 
row in the right-middle panel). [In the present case, 
Apg/ Ep = 0.012.] We have further verified that the size 
of the abrupt jump Z(~ 0.89) of the weight at kp in 
the repulsive case coincides (within 5%) with the value 
obtained from Fermi liquid theory [27| . which is related 
to the (retarded) self-energy according to the expression 
Z-i = 1 - [9Rc{S(fc,a.)}/aH,^,^_^^o; 

(iii) The large differences between the widths of the two 
structures of A{k,uj) in the repulsive case which persist 
across kp (left-lower panel) strongly deviate from the at- 
tractive case (right- lower panel), where the widths of the 
two branches reflect into each other at kp . 

It should be remarked that the coupling (kpap)^^ ~ 
— 1.5 wc have used to obtain the right panels of Fig. [11] 
is somewhat extreme, because very close to kp the sum 
of the widths of the two peaks of A{k,uj) exceeds their 
separation and the two peaks merge in a single one. To 
continue discerning two separate peaks in this narrow 
range of wave vectors near kp (specifically, from k ~ 
0.95fci? up to fc ~ l.5kp), a two-Lorentzian fit to the 
single broad peak is required. 

An example of this fit is shown at kp in the upper 
panel of Fig. [T71 where the peak at higher energy (short- 
dashed line) has larger weight than the peak at lower 
energy (long-dashed line) since kp > k^ (consistently 
with the right panels of Fig. [TB)) . Typical values of the 
X^-fit to isolate the two Lorentzians do not exceed 10~^. 
It is through this kind of fit that we were able to identify 
the value of the pseudo-gap Apg reported above, even in 
this rather extreme situation. 

Although in the attractive case near kp the two peaks 
of A(k,uj) merge apparently into a single one as a con- 
sequence of their broadening, the single-particle density 
of states N{Ld) given by Eq.® still maintains a well- 
pronounced feature about = due to the underlying 
pseudo-gap. This is evidenced by the dip occurring in the 
dashed line of the lower panel of Fig. [T7] corresponding to 
the attractive case, whose width is about O.li^i?- By con- 
trast, the dip is absent in the full line in the lower panel 
of Fig. [iTl which corresponds to the repulsive case and 
reproduces near = the frec-fermion result per spin 
component (dotted line) whereby A(k,uj) ~ S^lj — ^k)- 

The identification of the pseudo-gap is more direct 
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FIG. 17: Upper panel; Two-Lorentzian fit (short- and long- 
dashed lines) of the single broad feature of A{k = kF,ijj) (full 
line) in the attractive case with (/cfOf)"^ = —1.5 and T = Tc- 
Lower panel: Density of states vs uj for an attractive Fermi 
gas at T = Tc with (kpaF)^^ ~ —1.5 (dashed line), and for a 
repulsive Fermi gas at T = with (AfOf)"^ ~ 1-5 (full line). 
The values of vq and fco are the same of Fig. 1161 The dotted 
line shows the result for a free Fermi gas. 
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FIG. 18: Dispersions, weights, and widths at T = Tc ex- 
tracted from the two peaks of A{k,Lu) about kh/kp = 0.997± 
0.001 in the attractive case with (kpap)^^ ~ —0.8. Energies 
are in units of Ef- The inset shows the profiles of A{k,Ljj) 
vs Lj for three wave vectors k/kp = (0.980, 0.992, 1.004) that 
correspond to dashed, full, and dotted lines, in the order. 



when considering values less extreme than (kpap)'^ = 
— 1.5 for the attractive coupling, for which two distinct 
peaks in A{k,u!) can be distinguished even near k^. This 
is shown in Fig.[T8]ncar kp for the coupling {kpap)^^ = 
—0.8 still on the BCS side of the crossover, for which 
we get Apg/Ep — 0.127 ± 0.005. As already mentioned, 
even when the two peaks can be clearly distinguished in 
the region of the avoided crossing, a two-Lorentzian fit 
improves the accuracy of the determination of the disper- 
sions, weights, and widths associated with these peaks. 
Note from Fig. [TH] that, already at this less extreme cou- 
pling, moderate deviations from a simple BCS form arise 
in the dispersions and weights. 



IV. CONCLUDING REMARKS 

In this paper, we have focused a great deal of our atten- 
tion on the "pseudo-gap physics" which results near kp 
from the effects of attractive pairing fluctuations above 
Tc- From a careful analysis of the single-particle spectral 
function A{k,u}), we have identified the essential charac- 
teristics of A(k,uj) that can meaningfully be transposed 
above Tc starting from a mean-field description below 



Tc, even in situations when the inter-particle coupling is 
quite strong. Accordingly, the occurrence of a pseudo- 
gap has enabled us to carry over above Tc in an approxi- 
mate way concepts and results that are firmly established 
below Tc, where a truly long-range order extends over the 
entire system. 

We have further contrasted the occurrence of the 
pseudo-gap near kp, which produces a characteristic 
back-bending of the dispersion relation associated with 
the peak of A{k,uj) at negative to, with the eventual evo- 
lution of this peak for k ^ kp. Along these lines, we 
have regarded the pseudo-gap energy Apg as being associ- 
ated with pair correlations established in the system over 
intermediate distances of the order of the inter-particle 
spacing kp^, while the contact physics k ^ kp re- 
sults from pair correlations of shorter range. We have 
thus shown that attractive and repulsive pair correlations 
have similar effects for A; 3> A:_f but yield drastically dif- 
ferent results near kp. 

To appreciate intuitively on physical grounds the oc- 
currence of pseudo-gap phenomena established by attrac- 
tive pairing fluctuations above the critical temperature 
Tc of the superfluid transition, it may be useful to draw 
an analogy with the occurrence of damped spin waves 
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which are present in ferromagnetic materials above the 
Curie temperature Tc when a strict long-range order is 
absent. As it is often the case when dealing with phase 
transitions, examples from magnetic transitions may help 
one envisaging related phenomena occurring in different 
kinds of transitions. 

In the case of a ferromagnet, what is lost above Tc 
is the long-range order which organizes the spins over 
the whole sample. Yet, the spins can remain organized 
over moderate distances, so that spin-wave excitations 
may still propagate over these distances. The lack of 
full long-range order manifests itself in the damping of 
these "local" spin waves, which broadens their frequency 
spectrum. Experimentally, it has been long known since 
the work of Ref . (sij that spin waves can survive in mag- 
netic materials above the transition temperature (see 
also the more recent work of Ref.js^l and the refer- 
ences quoted therein). The accepted physical explana- 
tion for this effect is that a local magnetic order may 
exists above the magnetic transition temperature, such 
that spin waves can be supported over these length scales 
although with a shortened lifetime (this is especially true 
in low-dimensional systems - cf. Ref. [s^l)- Quite interest- 
ingly, experimental evidence has recently been collected 
that short-range spin waves may even underlie the mech- 
anism for high-temperature superconductivity [35 1. 

With this analogy in mind, pseudo-gap phenomena in 
a Fermi system with a mutual attractive inter-particle 
interaction can be envisaged as due to the persistence of 
a "local pairing order" above the superfluid temperature 
Tc, which takes place even when the (off-diagonal) long- 
range order is absent. The local order which is preserved 
by pairing fluctuations above Tc makes the single-particle 
excitation spectrum to resemble the one below Tc, al- 
though with an appreciable frequency broadening due 
to the temporal decay of these local excitations which 
cannot propagate over long distances in the absence of 
long-range order. 

In addition, this analogy may help one appreciating 
the connection between the presence of a pseudo-gap at 
low energy and the contact C at high energy [s^l, since 
they both derive from the existence of spatial correlations 
between fermions with opposite spins over moderate or 
else short distances. 

It is, finally, relevant to recall that the persistence of 
spin waves in the normal phase of magnetic materials had 
spurred a strong debate in the literature, their evidence 
being disputed in favor of a broadening of the energy dis- 
tribution [3^. Nevertheless, the existence of spin waves 
in the normal phase was eventually accepted by the oc- 
currence of new direct experimental evidence [33| . 

Something similar is apparently going on at present for 
the occurrence of a pseudo-gap in the normal phase of a 
Fermi gas with attractive inter-particle interaction, its 
evidence from momentum-resolved photoemission exper- 
iments P, 0, Q having been disputed in favor of a more 



conventional Fermi-liquid picture from thermodynamic 
measurements [37| . Similarly to what happened for spin 
waves above Tc in magnetic materials, the current dis- 
pute about pseudo-gap phenomena above Tc could possi- 
bly be brought to an end by a new evidence coming from 
direct detection of the upper branch of the dispersion of 
A{k,u)). This should be possible for a Fermi gas with 
attractive interaction via the experimental technique of 
Ref. [H once the temperature is raised sufficiently to pop- 
ulate the upper branch or for condensed-matter sys- 
tems where hole-like bands may give similar access to 
unoccupied state. 

This work was partially supported by the Italian MIUR 
under Contract Cofin-2007 "Ultracold Atoms and Novel 
Quantum Phases" . 



APPENDIX: T-MATRIX APPROXIMATION FOR 
A DILUTE FERMI GAS WITH REPULSIVE 
INTERACTION 

For a repulsive interaction of strength i;o > 0, the cut- 
off fcp in Eq.(IT]) cannot be let — )■ 00 in order to have a 
(positive) non- vanishing scattering length ap. For this 
reason, a potential of finite range (« k^^) has to be re- 
tained, without reaching the limit of a truly contact po- 
tential. In this case, one expects quite generally ap to be 
small, of the order of k^^ [38|. Solving for fco in Eq.([T]) 
one obtains: 



2ap 



2 7r" 

mvo 



(19) 



This relation sets the upper bound tt / (2kpap) on ko/kp, 
which is reached for very large values of mkpVQ. In prac- 
tice, we are interested in determining the behavior of the 
structures of the single-particle spectral function A{k,u}) 
for k fa kp or larger, in such a way that ko should ex- 
ceed these values of k. On the other hand, too large 
values oimkpvo (exceeding, say, 20) make the numerical 
summation over Matsubara frequency in Eq.(IS]) exceed- 
ingly difficult. As a consequence, a compromise has to be 
reached about the values of the coupling (kpap)~^ and 
mkpVQ. 

A systematic way to deal with finite values of the cutoff 
/co and the positive interaction strength vq is to introduce 
a separable potential in fc-space, of the form: 



y(k,k') = 1-0^(^0 -|k|)0(A:o-|k'|) 



(20) 



9{k) being the Heaviside step function. In this case, the 
particle-particle ladder reads: 

ro(k, k'; q, n,) = e{ko - |k|) e{ko - |k'|) fo(q, n,) (21) 

where 



1 



ro(q,i}.)-i = - 



dk 



«o J (2^)3 



0(fco - |k|) 



(22) 
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takes place of the expression (|6|) that holds in the limit 
ko — >■ CO. 

Finite values of ko and vq affect also the asymptotic 
behavior of ro(q, for large |q| and which is now 
given by: 



fo(q,f^.) 



«(q) 



(23) 



with 
a(q 



6{ko - |k|) [1 - /(a) - /(Cq-k)] 



(24) 



° I (27r)3 
The self-energy 

S(k,c^„) = -0(fco- |k|) y ^fcsr^ fo(q,17,) 
X e'"-" Go(q - k, - u;„) [2 - e{ko - |q - k|)l25) 

(with the distinct Hartree and Fock contributions) can 
then be conveniently organized as follows: 

E(k,w„) = 9{ko - |k|) [So + S'(k,w„)] . (26) 

Here, Eg = v^Uf /2 with Uf given by Eg. pop in the zero- 
temperature limit of interest and for fcp > ■\/2m/i. while 

E'(k,c^„) = - y^fcsT^ [fo(q,r!,) 



Go(q - k, r!, - uJn) [2 - 0(fco - |q - k| 



(27) 



where the convergence factor e^^""^ has been dropped ow- 
ing to the ri^-decay of the summand for large 

We are once more interested in the limit when k^/(2m) 
or \uJn\ are much larger than fc^T and fi. In this limit, 
the self-energy p6| reduces to the form: 



S(k,t^„) 



■9{ko^\M) 



-n^to(k,^„) + ALGo(k,- 



(28) 



in analogy to Eq.([9]), where is still given by Eq.(fT2|) 
with fo(q, fl^) replacing ro(q, ftiy). 

To obtain the single-particle spectral function accord- 
ing to Eq. ([7]) , analytic continuation za;„ — > w -I- zry to the 
real frequency axis needs to be performed. The imagi- 
nary part of the retarded self-energy is now given by: 



ImE(k,a;) 



^(fco-|k|) 



(2^ 



[2 - 0{ko - |q| 



X Imfo(k + q,a. + eq) [/(Cq) + + Cq) 



(29) 



where b{E) = [c^/C^bT) _ i ^j^g Bose function, in 
terms of which the real part ReE(k,w) can be obtained 



via Kramers-Kronig transform. To obtain the above ex- 
pression we have used the spectral representation 



fo(q,f^.) 



+ 00 



duj Imfo(q, w) 



(30) 



which generalizes to the present context an analogous 
expression considered in Ref.Q for an attractive zero- 
range potential. 

It is interesting to note in this context that /(^q) + 
6(^q) = in the zero-temperature limit, so that from 
Ea. p9)) one obtains ImE(k, = 0) = in that limit. 
The vanishing of ImE(k,aj ~ 0) is characteristic of a 
Fermi liquid at zero temperature for which a repulsive 
interaction is appropriate, and distinguishes it from a 
Fermi system with attractive interaction where superflu- 
idity sets in at the critical temperature Tc. A comparison 
of ImE(k, w) over an extended range of w about w = 0, 
for a repulsive interaction at T = and an attractive 
interaction at Tc, is shown in Fig. 1191 when fc = |k| = kp 
for the characteristic couplings {kpaF)~^ — ±1.5, in the 
order. 



-0.02 




FIG. 19: Imaginary part of the retarded self-energy vs uj at 
k — kp for a Fermi gas with repulsive (upper panel) and 
attractive (lower panel) interaction. Upper panel: T — 0, 
{kpap)'^ = 1.5, rnkpyo = 20, and kg — l.STkp- Lower 
panel: T ^ Tc = 0.0576rF and (kpapy^ = -1.5. 



These differences in ImE(k,CL>) between the repulsive 
and attractive cases result in marked differences in the 
single-particle spectral function A(k,uj), as shown in 
Fig. [201 over an even more extended range of uj for three 
distinct values of k about kp and with the same param- 
eters of Fig. [H 

A comment about the structure of A(k, uj) at nega- 
tive UJ for large k is in order at this point. Quite gener- 
ally, ImFo (k 4- q, a; + ^q) in the expression is non- 
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FIG. 20: Imaginary part of the retarded self-energy (left pan- 
els) and single-particle spectral function (right panels) vs u 
for three values of k, corresponding to a Fermi gas with repul- 
sive (full lines) and attractive (dashed lines) interaction with 
the same parameters of Fig. 1191 The log scale in the right 
panels highlights the small weights of the tails. 



vanishing provided the condition 

(k + q)2 



4m 



2/i 



(31) 



is satisfied. In addition, when a; < the sum of the 
Fermi and Bose functions in Eq. ip^ requires that yu < 

2 

^ < /i — w. When uj < 0, both conditions imply that 

^^^^3^ 2fi < 0. For large |k|, this means that |q| should 

be comparable with |k|, and therefore that ImE(k,w) 
of Eg. ([29|) is non- vanishing only for ~ — /i. A 
more stringent condition on the frequency interval where 
ImS(k,w) 7^ (and thus ^(k,w) 7^ 0) for w < and 
large |k| was considered in Ref. jll|. 

Some additional technical differences, between our ap- 
proach and the treatment of Ref. [Ill on the large-fc be- 
havior of the single-particle spectral function for a Fermi 
gas with repulsive interaction, need also be mentioned. 
In that work, emphasis was given to the large-A: be- 
havior of A{k,uj) in the repulsive case, to show that a 
structure at negative uj develops also in this case whose 
spectral weight is proportional to the contact C. For 
that purpose, it was sufhcient to consider coupling values 
(kpap)^^ ^ +10 corresponding to an extremely diluted 
situation, for which the details of the inter-particle po- 
tential of short range are irrelevant. In that limit, it was 
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FIG. 21: Single-particle spectral function of a Fermi gas with 
repulsive interaction {kpap)'^ = 4.4 at T = for k = A.Skp 
(here, ko = bkp and rnkpyo = 10). The broad peaks from 
the complete calculations (full lines) are contrasted with the 
delta-like spikes (dashed lines) resulting when in the self- 
energy the large-fc limit is taken before analytic continuation 
to real frequencies. 



possible in Ref.[ll| to use for the particle-particle ladder 
an expansion of the expression ([5]) up to second order in 
ap that would hold, in principle, only in the attractive 
case. For these rather extreme coupling values, the re- 



sults obtained in Ref.[ll| coincide with ours, which are 
obtained through a separable potential with parameters 
ko and vq related via Eg. ([TO)) . We recall that our use 
of a separable potential is required by the fact that we 
have extended the calculation of A(k,u;) in the repulsive 
case down to much smaller values (« +1) of the cou- 
pling {kpaF)~^ , in order to be able to compare with the 
corresponding results in the attractive case when k kp- 

As a final comment, we warn that, by performing the 
analytic continuation to the real frequency axis directly 
on the asymptotic form ([25)1 of the self-energy for large 
k < ko, one would get an expression of the type (jl3p for 
the single-particle spectral function with a sharp peak 
at w = — ^k- The correct procedure, of first performing 
the analytic continuation on the self-energy ([27|) and only 
then taking its limit for large fc, results instead in a broad 
structure of the single-particle spectral function at nega- 
tive frequencies, although with the same area A^/(4^^) 
of the approximate form ([T5|) for the attractive case. 
This is explicitly shown in Fig. [21] for k = A.5kp when 
(kpap)^^ = 4.4. This kind of non-commutativity, be- 
tween taking the analytic continuation and performing 
the large- |k| expansion on the self-energy, was already 
pointed out in Ref.ji^l while extending the t- matrix ap- 
proach to the superfluid phase below Tc- 

Although for large k the feature in A{k,cj) at neg- 
ative u! has formally the same weight A^/(4^^) in the 
attractive and repulsive cases, the numerical value of Aco 
(and thus of the contact C) depends on the sign of the 
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interaction and differs in the two cases. For instance, 
we obtain Aoo/Ep = 0.0860 for the attractive case and 
^oo/Ep = 0.1102 for the repulsive case using the pa- 
rameters of Fig. [121 These values can be compared with 
those obtained analytically from the Galitskii theory of 
a dilute Fermi gas [s^, according to which to leading 
orders in kpap'- 



A, 



27rn|aF| 



1 



(11 -2\n2)kFaF 



(32) 



where only the second term within braces depends on the 
sign of ap. This expression yields Aoo/ Ep = 0.0835 for 
the attractive case and Aoo/ Ep ~ 0.1063 for the repul- 
sive case, values which compare well with those obtained 
above from our numerical calculations. 
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